Anti-ferromagnetism, spin-phonon interaction and the local-density approximation in 

high-TV superconductors. 
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Results from different sets of band calculations for un- 
doped and doped HgBa2Cu04 show that small changes in 
localization can lead to very different ground states. The nor- 
mal LDA results are compared with 'modified' LDA results, 
in which different linearization energies make the O-p band 
more localized. The ground states in the normal calculations 
are far from the anti-ferromagnetic ones, while nearly AFM 
states are found in the modified calculations. The proxim- 
ity of an AFM state in the doped system leads to increased 
X s f, and the modified band structure has favorable conditions 
for spin-phonon coupling and superconductivity mediated by 
spin fluctuations. 



The Fermi surfaces (FS) of high-Tc superconductors 
contain barrel-like cylinders, which originate from the 2- 
dimensional CuO-planes. Other smaller FS-pockets or 
ridges are found in the more complex high-Tc materials, 
but it is probable that superconductivity is localized to 
the barrels, since they are common to all of the high-Tc 
cuprates. Electronic structure calculations, based on the 
local density approximation, LDA 1 , describe the FS cor- 
rectly in the doped cuprates. But the anti- ferromagnetic 
(AFM), insulating ground state of many undoped mate- 
rials is absent 2 . Thus, the band structure of the doped, 
metallic and superconducting system can be described 
within LDA, while the competing AFM insulating state 
in the undoped case is too high in energy to enter into 
the picture of possible ground states. This is a major 
problem when using LDA (and gradient corrected den- 
sity functionals, GGA 3 ) for many oxides, and the con- 
sequence is that most theories for high-Tc superconduc- 
tivity are oriented towards models of the band structure. 
This is generally done through introduction of parame- 
ters in so-called t- J-models, or through models for strong 
on-site correlation. 

It is not satisfactory to abandon the density functional 
(DF) theory without knowing the reason for this partial 
failure, especially since it is working so well for the large 
majority of other materials. Correlation is often cited as 
a possible reason, but why should this suddenly become 
large because of a slight change in chemical potential via 
different doping? The energies involved in on-site corre- 
lation, U, are large, several volts are typical, so U has to 
vary extremely fast with doping. Here I suggest another 
possibility; that rather small changes of the potential can 



make that two very different configurations have similar 
total energies and enter into competition. Thus, the sit- 
uation could be like the one for the ground state of iron, 
where LDA calculations predict nonmagnetic fee as the 
ground state, while GGA correctly predicts the FM bec 
state as the stable one 4 . The differences in total energy 
between the different configurations are small within each 
DF-potential, only the fine details make one state more 
stable than the other. Still, the changes in potential given 
by GGA are small, and many other properties of Fe can 
be calculated correctly by the use of both DF-potentials. 
The GGA applied to oxides make no large improvement 
compared to LDA. 

Band calculations, using the linear Muffin-Tin Orbital 
method (LMTO) 5,6 and LDA for one of the high-T c ma- 
terials, HgBa2Cu04 (HBCO), showed that spin fluctu- 
ations are coupled to phonons within the CuO plane 8 . 
The qualitative results for the doped case are consistent 
with many observations in high-Tc materials such as the 
existence of stripes, pseudogaps in the density-of-states 
(DOS), isotope shifts, AFM fluctuations, phonon soften- 
ing in the doped materials etc.. But the amplitude of the 
coupling is too low for a high Tc- The likely problem 
with LDA is that the total energy of the AFM ground 
state is much above the nonmagnetic state in undoped 
HBCO. We shall analyze the ingredients of the total en- 
ergy of an AFM configuration by comparing it to the case 
with no spin-polarization. The exchange energy is always 
in favor of spin-polarization, as is well known from the 
Stoner model for FM 9 . In contrast to the FM case, also 
the kinetic energy favors magnetism if the AFM state 
will open a gap or pseudogap at the Fermi energy, Ep, 
as was discussed previously for HBCO 7 . But the AFM 
state, where the potential for one spin show a spacial 
modulation, costs hybridization energy. The tail of a 
wavefunction of one spin from one site enters another site 
where the spin is different, making the conditions for hy- 
bridization different from the unpolarizcd case. The FM 
case is well described by a rigid band shift of the entire 
majority band structure relative to the minority one, but 
in the AFM case there is no such shift. One can expect 
that the hybridization would decrease for more localized 
states, since the tails become smaller. The plane oxygens 
(between Cu sites of opposite spins) have their p-states 
hybridizing with the d-states on Cu, and vice-versa. The 
idea is to make the O-p states more localized so the in- 
teraction between Cu of opposite spins can decrease and 
thereby minimize the cost of hybridization energy for an 
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AFM state. 

The logarithmic derivatives, Dg, of the LMTO 
method 5 , provide information about localization of each 
band I. The values of Dg are negative in the band region, 
and they decrease continuously for increasing energy un- 
til they become infinite at the top of the ^-band. The 
localization is best when the energy of the state is at the 
top of the band. The LDA-value of Dg for the oxygen p- 
band in HBCO is about -5 at Ep , showing that the band 
is not completely filled. The hybridization energy should 
diminish when D t — > — oo, corresponding to a very local- 
ized state at Ep. States far below Ep are extended, but 
such states are less affected by the spin splitted poten- 
tial. In figure 1 it is seen that the nonmagnetic and AFM 
DOS functions differ only in the upper part of the Cu-d 
band. This is because the radial extent of the exchange 
splitting of the potential is localized within the Cu-atom, 
and the splitting of an extended state is relatively small. 

The optimal choice of the linearization energies for the 
i'-bands in LMTO is at the center of each band. One way 
of making the states more localized is to do the lineariza- 
tion at lower energies, at the bottom of the bands. This 
makes the self-consistent -D^-value at Ep positive, show- 
ing that the top of the O p-band goes down compared to 
the normal case. This is the main difference for making 
modified (more localized) bands without changing the 
LDA-potential. It is not a 'good' way of calculating the 
LDA bandstructure, but as will be shown, the results 
serve as examples of how to generate an AFM, insulating 
state in the undoped case and stronger AFM fluctuations 
in the doped case. A full gap appears in the DOS of the 
undoped HBCO if a staggered field of ±23mRy is ap- 
plied within the Cu sites in the normal LDA calculations. 
Calculations with lower choices of linearization energies 
('modificd-I') give gaps already at smaller applied fields, 
~10 mRy, see figure 1. Stronger changes are found when 
an additional modification fmodified-H') is made in the 
calculation of the Madelung potential. The method for 
calculating the Madelung potential uses non-overlaping 
MT-spheres in the calculation of the point charges, that 
are smaller than the Wigner-Seitz (WS) spheres used in 
the LMTO scattering matrix. An interstitial volume be- 
tween the MT and WS spheres has an interstitial charge. 
This procedure has a moderating effect on the Madelung 
potential for compounds 6 . An additional lowering of the 
p-band is observed when using oxygen MT- and WS- 
spheres of equal size. For this 'modified-IF LDA calcula- 
tion it is just possible to obtain a gap and AFM without 
an applied field, sec figure 2. The moment is 0.23 /xg per 
Cu atom. 

The band widths of the d- and p-bands from the CuO- 
plane are reduced by 0.02 Ry in the 'modificd-IF LDA 
calculation. A more important change is that the top 
of the O-p band goes down from 0.14 Ry above Ep to 
about 0.17 Ry below Ep. The top of the Cu-d band goes 
up relative to Ep, by 0.03 Ry. This can be seen in the 
DOS, in which the shoulder of large DOS moves upwards 
closer to Ep. The same happens when large fields are 



applied and the shoulder can even reach the position of 
Ep for the doped case, cf figure 1. This is an extreme, 
not realistic case and it will hardly help the AFM state. 

The charge distribution is changed in the modified 
case. As can be expected from the changes in band 
positions there is an increase of the the O p-charge by 
0.11 e~, while the Cu d-charge decreases by 0.09 e~ 
per atom. These changes follow the same trends that 
were found from calculations with improved correlation 
added to the LDA for the CuO planes in the 'infinite- 
layer' structure 10 . No detailed comparison with results 
from the method of Stollhoff is possible because of differ- 
ent basis, different structure, and so on. But the largest 
charge redistribution caused by correlation 10 , a decrease 
of the x 2 — y 2 Cu d-charge by 0.38 e~ and an increase 
of the O p-charge along the bonding by 0.19 e~, are in 
line with the present results. The result of Stollhoff call 
for moderate corrections of LDA due to correlation, but 
it also excludes the limit of strong correlation and the 
Mott-Hubbard model 10 . 

Before turning to the doped system it might be appro- 
priate to make some comments about the energy scale 
of the perturbations. The energy of the magnetic field, 
/iH=23 mRy, required to make the undoped system an 
AFM insulator in LDA calculations, seems large in view 
of the fact that 1 mRy corresponds to a magnetic field 
of 230 T. However, this is not much on the scale of other 
electronic energies in the band calculation. Total ener- 
gies calculated in LDA and GGA can differ by a few Ry 
per atom, although differences in total energies between 
different configurations (that are to be compared with 
physical properties) and other main features of the band 
results are reasonable from both DF potentials. Thus, 
the value 23 mRy can by itself not be taken as an indica- 
tion that LDA is very far from a good description of the 
ground state. The fact that fxH can approach mRy and 
produce an AFM insulating state in 'modified' LDA cal- 
culations, shows that electronic energies associated with 
charge localization can easily replace the energy coming 
from an applied magnetic field. Moreover, it can be re- 
called that when LDA bands are replaced by band models 
including strong correlation one needs to apply larger en- 
ergies in terms of on-site correlation U, which typically 
are of the order of 5 eV (~400 mRy). 

The next step is to do the same modification for the 
hole doped case. The pseudogap will appear at the new 
position of Ep if the AFM state is modulated, as has 
been shown from calculations using the virtual crystal ap- 
proximation for a supercell oriented along the CuO bond 
direction 7 . The wavelengths of spinwaves and phonons 
differ by a factor of two for an optimal spin-phonon 
coupling, and they depend on the doping. This, and 
other features of the spin-phonon coupling can be found 
in previous publications 7 ' 8 , but one may note that the 
relatively short supercell considered here corresponds to 
phonon and spin waves in an overdoped system. The 
pseudogap, the large exchange enhancements and mag- 
netic moments on the Cu-sites, are all indications of a 
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metal-insulator transition. The normal LDA exchange 
enhancement, S, evaluated as the exchange splitting of 
Cu-d band divided by the applied field, is about 2.0 when 
the cell contains the phonon distortion and the spinwave. 
The enhancements in the 'modified' LDA calculation are 
much larger, see Table 1. The moments (m) are larger 
than in the normal LDA, and they are normally propor- 
tional to the applied field, fiH, to = M ■ fj,H. 

The calculation of A s / is made assuming harmonic con- 
ditions for a 'frozen' spinwave 8 ' 11 : 

A s/ = NI 2 /2F (1) 

where N is the paramagnetic DOS at F,p (100-115 
states/cell). The matrix element, /, is the derivative of 
the potential with respect to m. In the harmonic limit; 
< SV >= I ■ m, and / is calculated from the FS average 
of the difference of band energies, e™ — e , between the 
cases with and without magnetization. It is assumed that 
the difference in total energy between the configuration 
with moment m and at zero moment, E m — Eq, depends 
almost quadratically on m, E m — E = Fm 2 , similar to 
the dependence of the total energy on the atomic dis- 
placements, Ar, for harmonic phonons. 

Table I shows the results of N, S, M, I 2 , F and A s / 
compared to the results from the ordinary LDA calcula- 
tions, case (b) in the previous work 8 . The enhancement 
S and the factor M show qualitatively how the system 
gets closer to an AFM transition. The calculations are 
made for applied fields up to 50 mRy in the ordinary 
LDA, while up to 30 mRy in the modified cases. The 
moments are larger in the modified calculations, when 
the AFM instability is closer and the harmonic behav- 
ior is less evident. The proximity of the AFM instability 
makes the ingredients (I 2 and F) to A s / smaller, and the 
relative errors in F and A s / are larger. It is difficult to 
calculate A s / in the limit m — > 0, since an estimation 
of F from the differences of total energies between small 
and zero fields is too uncertain. Using the high or low 
TO-limit for X s f makes no big difference in the ordinary 
LDA, but in the modified calculations it is the low-m 
limit that gives the largest A s /. The results displayed in 
Table I are based on results from larger to. Despite the 
increase in relative errors it seems clear that the coupling 
constant for spin-fluctuations is larger in the modified 
LDA calculations, and it can be large enough for a high 
Tc- However, the calculation of Tq requires further con- 
siderations, whether the prefactor in a BCS-like formula 
should be given by the excitation energy from phonons 
or from spin waves, whether the corrections due to non- 
harmonicity are important, etc.. No theory is yet known 
to describe superconductivity within a spin-phonon cou- 
pled system, but it is reasonable to assume that a larger 
A s / should give a larger Tc- 

The main effect from the modified band calculation 
is that the potential becomes more attractive within the 
oxygens compared to the normal LDA. The HBCO struc- 
ture contains Ba, Hg and apical oxygens, and the influ- 
ence from these sites on the potential within the CuO 



plane gives additional possibilities for AFM. For instance, 
it was found that the relative positions between the CuO 
plane and Ba or apical oxygens are important for the 
local exchange enhancement on Cu 7,8 . It remains to be 
investigated in detail how real phonons with longer wave- 
lengths interact with spin waves in underdoped systems. 
Various reasons can be suggested to lead to a more at- 
tractive potential on oxygen sites, such as non-locality, 
kinetic energy or correlation, but this will not be dis- 
cussed here. 

In conclusion, two band calculations using the same 
LDA potential, but different linearizations, show very 
different results for the stability of the AFM state of un- 
doped HBCO compared to ordinary LDA. The same ex- 
ercise for doped HBCO show that the modified LDA will 
enhance the spin-phonon coupling and A s /. It is argued 
that LDA might need rather modest corrections in or- 
der to promote an AFM, insulating state in oxides. The 
resulting charge transfer is in agreement with an indepen- 
dent investigation of moderate correlation 10 . The gap is 
in this case an ordinary band gap, caused by modulations 
of AFM within the CuO plane. This is in contrast to the 
Mott-Hubbard gap due to very strong correlation. 



1 W. Kohn and L.J. Sham, Phys. Rev. 140, A1133, (1965); 
O. Gunnarsson and B.I Lundquist, Phys. Rev. B13, 4274, 
(1976). 

2 T. Jarlborg and G. Santi, Physica C 329, 243 (2000), and 
references therein. 

3 JR. Perdew and Y. Wang, Phys. Rev. B33, 8800, (1986). 

4 J. Kubler, Phys. Lett. A81, 81 (1981); C.S. Wang, B.M. 
Klein and H. Krakauer, Phys. Rev. Lett 54, 1852, (1985); 
B. Barbiellini, E.G. Moroni and T. Jarlborg, J. Phys.: Con- 
dens. Matter 2, 7597 (1990), and references therein. 

5 O.K. Andersen, Phys. Rev. B12, 3060 (1975). 

6 T. Jarlborg and G. Arbman, J. Phys. F7, 1635 (1977). 

7 T. Jarlborg, Phys. Rev. B64, 060507(R) (2001), Phys. Lett 
A295, 154 (2002). 

8 T. Jarlborg, Phys. Rev. B68, 172501 (2003). 

9 T. Jarlborg and A.J. Freeman, Phys. Rev. B23, 3577 
(1981). 

10 G. Stollhoff, cond-matt/0211230 (2002). 

11 T. Jarlborg, Physica C 385, 513, (2003); Phys. Lett. A300, 
518 (2002). 



3 



TABLE I. Calculated parameters (see text) for ordi- 
nary and modified LDA for a spinwave co-existing with 
a halfbreathing phonon as in ref. 8 . The units are for N 
states/cell/Ry, for M: ^ B /mRy/cell, for I 2 : 10~ 3 Ry 2 / n 2 B , 
and for F: mRy / 'cell / ' \j? b . 
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FIG. 1. Calculated DOS of Hg2Ba 4 Cu 2 08 in which the 
linearization energies are low compared to the band centers 
('modified-I', see text). A staggered field of ±10mRy on 
the Cu makes a gap at Ef- Note that the AFM and the 
non-polarized DOS differ only at the upper part of the Cu-d 
band. 
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FIG. 2. Paramagnetic (bold line) and AFM (thin line) DOS 
near Ef of Hg2Ba4Cu20s calculated for the lowest choice of 
linearization energies and a modified Madelung charge ('mod- 
ified-IF, see text). 
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